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RECENT ADVANCES IN AERODYNAMIC ENERGY CONCEPT FOR FLUTTER SUPPRESSION 


AND GUST ALLEVIATION USING ACTIVE CONTROLS 
E. Nissim* 

Langley Research Center 


SUMMARY 

Control laws are derived, by using realizable transfer functions, which 
permit relaxation of the stability requirements of the aerodynamic energy con- 
cept. The resulting aerodynamic eigenvalues indicate that both the trailing- 
edge and the leading-edge — trailing-edge control systems can be made more 
effective. These control laws permit the introduction of aerodynamic damping 
and stiffness terms in accordance with the requirements of any specific system. 
Flutter suppression and gust alleviation problems can now be treated by either 
a trailing-edge control system or by a leading-edge — trailing-edge control sys- 
tem by using the aerodynamic energy concept. Results are applicable to a wide 
class of aircraft operating at subsonic Mach numbers. 


INTRODUCTION 

Flutter suppression and gust-alleviation systems using active controls 
tend to be very sensitive to system changes caused by different flight condi- 
tions (flight speed, flight altitude, flight duration, and type of mission). 

The aerodynamic energy concept (ref. 1) was formulated in an attempt to define 
•active control systems which do not exhibit such sensitivities to changing 
flight conditions. Although recent applications of the aerodynamic energy 
concept to specific problems of flutter suppression (refs. 2 and 3) and gust 
alleviation (ref. 3) yielded encouraging results, it was indicated that the 
derived control laws could be improved. 

In the original development of the aerodynamic energy concept, idealized 
transfer functions were employed. These idealized transfer functions were 
difficult to realize and only approximate implementations could be made 
(ref. 2). In addition, the aerodynamic energy concept contained stringent 
requirements of a sufficient, but not necessary, condition for stability which 
effectively ruled out the trailing-edge (T.E.) control as a single stabilizing 
control system in favor of a combined leading-edge — trailing-edge control 
system. 

The present work reviews the aerodynamic energy approach in light of the 
experience obtained since its original development and describes a way of relax- 
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ing the stringent stability conditions and improves the derived control laws by 
using realizable transfer functions. Both trailing-edge and leading-edge — 
trailing-edge control systems are investigated. 
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SYMBOLS 

leading-edge free parameter 
trailing-edge free parameter 
semichord length 

element i, j of control law matrix [c] 
element i, j of control law matrix [g] 
structural damping coefficient in ith mode 
bending displacement, positive in down direction 

= nR 

0)b 

reduced frequency, — 

Mach number 

mass ratio, defined in equation (B9) 
mass 

number of degrees of freedom 
energy dissipated by system per cycle 

generalized aerodynamic force along h and a, respectively 

transfer function, defined in equation (27) 

number of activated control surfaces 

wing semispan or reference length 

time 

flight speed 

oscillatory angle of attack of wing, positive in nose-up direction 

leading-edge and trailing-edge control surface deflections, 
respectively, positive in down direction 
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damping coefficient 
eigenvalue of [u] 


= k^X 


equivalent X for structural damping 
fluid density 


oscillatory frequency 


Matrices: 

[Ar] . [Al] 
[ar] , pi] 
[B] 

[c] 

[E] 

{F} 

Pr] » Pl] 

{Qo} 

{qR}.{qi} 

[T] 

P] 

[X] 

{Cr}» {C i} 


real and imaginary parts of aerodynamic matrix [aJ (see eq. (1)) 
real and imaginary parts of aerodynamic matrix [a] (see eq. (15)) 


inertia matrix 


control law matrix 


structural stiffness matrix 


column matrix of forces 


control law matrix 


real and imaginary parts, respectively, of energy eigenvector 
modal matrix 

complex response vectors 

complex amplitude of response vector 

real and imaginary parts of {(^q] 

transfer function matrix 

energy matrix 

diagonal matrix of eigenvalues of p] 

real and imaginary parts, respectively, of generalized energy 
coordinates 


eigenvector of p] 
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Subscripts: 


c 

control 


I 

imaginary part of a 

complex value 

max 

maximum 


min 

minimum 


R 

reference, also real 

part of a complex value 

s 

structural 


Other 

notations: 


1 1 

absolute value 


[]T 

transposed matrix 



complex conjugate 


LJ 

row matrix 


{ } 

column matrix 



Dots over symbols denote derivatives with respect to time. 

THE AERODYNAMIC ENERGY APPROACH 
Basic Concept 

An active control system on a lifting surface, such as a wing, actuates 
a control surface in response to oscillations of the wing in a manner which 
stabilizes the system. Because flutter instabilities originate from aerody- 
namic forces acting on the wing, the possibility arises of changing the aero- 
dynamic forces through an appropriate activation of control surfaces. The 
aerodynamic energy concept was developed (ref. 1) to examine this possibility. 
The concept centers around the work P done by aerodynamic forces on their 
surroundings, per cycle of oscillation. A control surface activation is sought 
that leads to large and positive values of P, thereby energy dissipation and 
stability are insured. 

Development of energy concept.- In this section, the energy concept is 
redeveloped (with details presented in appendix A) both to help the reader fol- 
low the present work and to show that the original derivation pertains only to 
mass-balanced control systems. 
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The n equations 


{f} = -0)2 [b + TTpb^s(AR + iAj)] {q} + [E]{q} 


( 1 ) 


represent 

controls, 

the equations of motion of 
where at flutter 

n structural modes with 

r activated 


{f} = 0 



and 




0) 

frequency of oscillation 



[B] 

n X (n + r) mass matrix (n 

natural modes and r active controls) 

[Ar] , [Ai] 

real and imaginary parts of aerodynamic matrix. 

respectively 

[E] 

stiffness matrix 



P 

fluid density 



S 

reference length 



b 

reference semichord length 



{q} 

response vector comprising 

n structural modes and 

r control 


modes, {q} = < > 

Iqc i 




The matrices in equation (1) can be partitioned into square matrices (n x n) 
relating to the structural modes (subscripted by s) and rectangular matrices 
(n X r) relating to control surface couplings (subscripted by c). 

After partitioning the matrices, equation (1) becomes 


{f} = ^-0)2 J^[Bsj bJ + TTpb^s^[AR^g| Ar^q] + i [A j , g; Aj ^ c])J + (2) 


Assume a control law of the form 


{qj = [T]{q} 


(3) 


where is a r x n matrix representing the transfer functions of the con- 

trol law. Substituting equation (3) into equation (2) yields 
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{f} = (-0,2|-|b3] . [bJ[t] . ^Pb^s([AR,3] . [AR,e][T] 
+ i[Al,s] + i[Al,c]M)] + [Eg] + [Ec][T]]{q} 


The matrix [Eq] is considered to be zero since it is assumed that no elastic 
couplings exist between the structural modes and the control deflections. As 
shown in appendix A, the work P done by the system on its surrounding per 
cycle can be written as (eq. (A12)) 


p = - 10lj|-([*I,s] • pl,j" • [Al,c]W * 


[Bj[T] - 

•jrpb^s 


^ [AR,s] - [Ar.s]" + [Ar,c][t] - [T»]’‘[Ar,c]^)|{qr 


where ( from eq . ( A6 ) ) 


{q} = {qo}e^“^ = {qR + iqile^^^t 


The sign of P determines stability, and therefore, it is advantageous 
to convert equation (5) to a more convenient form. As shown in appendix A 
(eq. (A20)), P can be reduced to the form 


p = -ir2pbVs([CRj[x]{CR} + LCiJ[x]{Cr 


P = - ^ + X2(cr,2 + ^ 1 , 2 ) + • 


/ 2 2 
+ ^nUR,n + Ci,n 


where ^Xj is a diagonal matrix of the eigenvalues Xj^, necessarily real, of 
the Hermitian matrix (eq. (Al4)) 
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M = -([Al.s] * * [Al,o] H * [T*]"[Al,o]’ 

* i([AR,s] - [Ar.s]" * [A«,o][T] - [T*]"[A„,o]" - 


MM - 

TTpb^S 


( 9 ) 


and where the vectors {?r} and {?i} are defined by the transformation 
(eq. (A15)) 

{Qo} = [Qr + (''0) 


The matrix [ISp + iQ£] is a square modal matrix of the principal eigenvectors. 

Discussion of e nergy concept .- The work per cycle P done by the sys- 
tem on its surroundings has a direct bearing on the stability of the system. 

If P is positive, the system is dissipative, and therefore stable. If P is 
negative, the system is unstable because work is done by the surroundings on 
the system. Equation (8) shows that if all the eigenvalues of the system 

are positive, the system is stable regardless of the motions represented by the 
generalized energy coordinates C . If one or more of the X eigenvalues is 
negative, the system is potentially unstable. Its ultimate stability is deter- 
mined by the relative values of the terms C and X, If the ^ values make 
the positive eigenvalues dominate the right-hand side of equation (8), the 
work P is positive and the system is stable. If, on the other hand, the 
^ values make the negative terms dominate equation (8), P is negative and 
the system is unstable. Hence, the requirement for all X*s to be positive 
is a sufficient but not a necessary condition for stability. 

For mass-balanced control surfaces ([BqJ = o), the eigenvalues X 
obtained from [u] (eq. (9)) are dependent only on the aerodynamic properties 
of the system and the activated control law ^matrix [t]). In the case of 
mass-balanced surfaces the eigenvalues are referred to as aerodynamic eigen- 
values and are, in general, functions of the reduced frequency k and Mach 
number M. If mass unbalance is a fixed quantity in the system, the eigen- 
values X also depend on the fluid density p in addition to their dependence 
on k and M. In the present work, only mass-balanced systems are treated and 
therefore eigenvalues are obtained from the following [[uj matrix (eq. (A21)) 

M = [-([Al,J * [AI,J" * [Al.JW ^ M^[Ai,J") 

a([ab,J - [ar.s]’' * [Ar.o][t] - [t*]’'[Ar,o]^)] 


Equation (11) shows how the transfer function matrix [t] affects the 
matrix [u] and therefore affects the aerodynamic eigenvalues. 
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Generalized Model 


The energy approach has been formulated for a general n degree-of- 
freedom system. Therefore, the energy concept can be applied to any problem. 
The results in this report, however, are specific for the system considered 
since the generalized aerodynamic forces depend not only on the system geometry 
but also on its structural modal responses. If, however, the energy concept 
is applied to a two-dimensional strip, the aerodynamic matrices are independent 
of geometry and responses of the system. As a result, the aerodynamic eigen- 
values are independent of any specific system and are only functions of k, 

M, and the transfer function matrix [t] . Therefore, if [t] is defined in 
accordance with the relaxed energy concept, using a two-dimensional strip as 
a model, these [t] values are applicable to any three-dimensional wing within 
the limitations of strip theory; thus, the model is generally applicable. In 
the present work, the two-dimensional strip model is adopted in much the same 
way as in the original derivation of the energy concept. Sketch (a) illus- 
trates the system considered, and the arrows indicate positive displacements 
and rotations. 


Undisturbed position 



Relaxation of Energy Concept 

The energy approach, in its original development, sought to determine 
the matrix [t] to render all the aerodynamic eigenvalues large and positive. 
This requirement regarding the aerodynamic eigenvalues insures both the sta- 
bility of the system (since P is always positive) and its insensitivity to 
various flight conditions (which manifest themselves in the form of changing 
values of X and changing values of the system responses C). Assume that 
a relaxation is now introduced which permits some of the aerodynamic eigen- 
values to be negative. Stability can be achieved under these conditions by 
modifying the responses of the system to render the responses associated with 
the positive eigenvalues to be the dominant ones. This latter requirement 
forms a necessary condition for stability but does not insure, in itself, the 
insensitivity of the resulting stabilized system to the various flight condi- 
tions. In order to Insure that this relaxed stability requirement yields a 
system which shows only small sensitivities to the changing flight conditions 
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the absolute values of the negative aerodynamic eigenvalues must always be made 
much smaller than those eigenvalues associated with the dominant responses of 
the stabilized system. For the generalized two-dimensional model adopted in 
this work, two aerodynamic eigenvalues, ^min ^max> obtained. In the 

original derivation of the aerodynamic energy concept, was required to 

be positive and large. In the relaxed energy approach, is permitted to 

be negative provided 


'^max 


l^minl 


(^min ^ 


and provided that this inequality is maintained for all flight conditions. 

This relaxation is made possible by abandoning the sufficiency condition for 
stability in the original formulation while maintaining its insensitivity to 
changes in flight conditions. It should be stressed at this stage that the 
generalized two-dimensional flutter model adopted herein serves only to indi- 
cate, on the basis of the strip theory, whether energy is dissipated or 
absorbed by the partial span strip where the activated controls are installed. 
Therefore, to suppress flutter with a minimum number of activated partial span 
strips, enough energy should be dissipated in the activated strip to compensate 
for energy input by the nonactivated portions of the wing. Not only should 
Xjnin be made positive but also Xjjjj^^ might assume large (and positive) values. 
Since the dissipation of energy by the activated strip depends both on Xjjjj^j^ 
and on X ^a y , the importance of X ^a y should not be overlooked even when 
Xmin is positive and large. Considerable improvements in the potential per- 
formance of the activated control system may result, if changes in the control 
gains are permitted which lead to small degradations in provided these 

degradations are accompanied by large increases in X ^g y . It is therefore 
required to determine the optimum values of the transfer function matrix [t] 
to lead to 


Xnjin = Near maximum value (may be negative) 


^max 


'•mini 


These two requirements regarding X^in ^max referred to in this 

work as the "relaxed energy requirements.” 

Formulation of a revised optimiz ation procedure .- The optimization proce- 
dure used in the original development of the energy concept consisted of the 
following basic steps: 

( 1 ) Define the generalized aerodynamic forces and acting on a 

unit span of the two-dimensional strip and express them in the form 
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irpbV 


([*R,s] * i[*I,s])[ ^ j* ([*B,o] * l[*I,o])|* j 


(12) 


where all the matrices are of order 2x2. 

r 3 ) r h/b) 

(2) Assume a form for the matrix [t] where ^ J - [t]^ j. 


(3) For a given Mach number M and initial values for the gains in [T], 
compute the numerical value of M for a specific value of k. 

(4) Compute the two eigenvalues resulting from [u] and designate them 

^min ^max" 

(5) For a fixed value of M and [t] , repeat steps (1) to (4) for differ- 
ent values of k within the range 0.0128 ^ k $ 19.5 (0.05 S 1/k S 78). 

(6) Determine the area under the curve defined by the variation of X mi n 
with 1/k. 

(7) Optimize the gains in matrix [t] to obtain a maximum area under the 
curve of X mi n a function of 1/k. 

This seven-step procedure, which used the area under the curves of X mi n 
against 1/k as a target function, was found to be generally satisfactory. 
However, a plot of X^j^n against 1/k and X ^ay against 1/k for the 
unactivated system (fig. 1) shows that the very low frequency part of the curve 
contributes most to the target function. It is therefore possible that large 
improvements in the values of X mi n at intermediate frequencies may be com- 
pletely "masked” during the optimization process if accompanied by very small 
degradations of ^min at very low frequencies (high 1/k values). 

A reexamination of the X „n~ n against 1/k behavior, undertaken in the 
present work, has led to the following two initial changes: 

(a) The very large aerodynamic eigenvalues at the very low k range fol- 
low from the representation of the aerodynamic forces 


{Fj^} = irpb^a)2s^[AiJ + i[Ai]^{q} 


(13) 


The matrices [Ap] and [Aj] do not contain the true variation of the aero- 
dynamic forces with k because of the terms left outside the matrices. 

The frequency oo is not constant and varies with every mode of oscillation. 

A more rational representation of {FaI > which maintains the full dependence 
of the aerodynamic matrices on k, can be expressed by 
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{Fa> = Trpb2v2s([AR] + i[5i]){q} 


(14) 



where V is the flight speed and 

[Ar] + ipi] = k2^[AR] + ipi]) (15) 

Equation (14) leads to the following change in equation (8) (for n = 2): 


P = - ■»r2pb2v2s[jCi(cR j + 5l,l) + ^2^T{,2 + ^ 1 , 2 ) (16) 

where 

Xi = k2\i (17) 

Hence, for k < 1, the newly defined eigenvalues are smaller than the 

originally defined eigenvalues Xj^ by a factor of k2. 

(b) The k range in the original derivation (0.0128 S k ^ 19.5) was 
unnecessarily wide. In the low k range, gust-alleviation problems generally 
occur around k » 1/8; in the high k range, the accuracy of the aerodynamic 
derivatives is questionable for values of k > 3.5. Therefore, the k range 
adopted was 

1 

— ^ k S 3.5 
25 

or 

1 1 

S — $ 25 

3.5 k 


The optimization algorithm used in the original development of the energy 
approach was found to be deficient in many respects. The algorithm used in the 
present work consists of a variation of Stewart adaptation (ref. 4) of the 
Davidon-Fletcher-Powell method (ref. 5). The variation introduced permits the 
different elements of the matrix [tJ to be constrained within preassigned 
limits without resorting to penalty functions and with excellent convergence 
characteristics . 

In accordance with the relaxed energy concept there remains the problem 
of permitting small degradations in if accompanied by large increases 

X may . The initial approach was based on adding to the target function a 
weighted value of the area under the curve of X „,aY against 1/k. This 
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approach, although effective, leaves the value of the weight as an arbitrary 
parameter for the designer. However, as the present work progressed, the 
parameters which appreciably increase ^ax with relatively small degradations 
in X mi n have been identified, and the optimization procedure was modified in 
a manner which avoided weighting X^ax* This point is discussed further in a 
later section of this paper. 

Data and scope of optimization.- As in reference 1 and indicated in 
sketch (b), the two-dimensional strip has constant 20-percent-chord leading- 
edge (L.E.) and trailing-edge (T.E.) control surfaces. Similarly, a reference 
point for sensing the motion of the main surface was kept constant at 30 per- 
cent of the chord (measured from the leading edge) . The gains associated with 
the control laws presented in this work refer therefore to the 30-percent chord 
point. These gains can be modified to accommodate a sensor located at a dif- 
ferent point along the chord by employing a simple transformation matrix 
(ref. 1). 


Undisturbed position 



The investigation covers both a trailing-edge and a leading-edge — 
trailing-edge control system. Each of these control systems is driven by two 
different realizable transfer functions represented by the matrix [t] . The 
gains associated with each transfer function are optimized within a predeter- 
mined range of values by following the relaxed energy approach. Generalized 
control gains were tested for adequacy over wide ranges of reduced frequen- 
cies k, of subsonic Mach numbers M, and of maximum permissible gain values. 


PRESENTATION AND DISCUSSION OF RESULTS 

Results for the unactivated strip at various subsonic Mach numbers are 
presented first. These results form the basis for assessing any improvements 
introduced by the activated control systems. The results relating to each of 
the two transfer functions employed in this work are presented and discussed 
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separately. For each transfer function, a study is made of both the trailing- 
edge and the leading-edge — trailing-edge control systems. 


Basic Nonactivated System 

The variation of X^ln ^max with 1/k at different values of M, 

for a wing strip having no control surfaces (that is, [T^ = o) is shown in 
figure 2. The quantity Xmin negative through most of the 1/k range, 
whereas Xmax positive throughout the range of 1/k. Furthermore, the 
negative Xmin values are of the _same order of magnitude as X^^x* Note also 
that the aerodynamic eigenvalues X presented in figure 2 are greatly reduced 
in the high 1/k range when compared with the original eigenvalues X 
(fig- 1). 


Damping Type Transfer Function 

The damping type transfer function (D.T.T.F.) is similar in form to the 
one used in reference 3; that is, 


or 


3 ' 

6 




h/b' 


a 


(18a) 


• 3 ' 

. 6 , 


([C] * 


(18b) 


where [c] and are 2x2 constant matrices. Therefore, [t] is given 

by 


[T] = [C] + Ik [G] 


(19) 


The control law represented by equations (18) differs from the one used in the 
original derivation of the energy approach ( which used the term i [g] instead 

b \ ' b 

of the — iw[G] term j. The implementation of the — iw[G]] term in equa- 
tion (l8a) is straightforward, whereas the implementation of the i[G] term 
could only be achieved with great difficulties and in an approximate manner. 

The damping type transfer function was chosen for the relaxed energy approach 
because it is similar to the original control law and has been shown to be 
effective when applied to specific configurations (refs. 3 and 6) ^by using the 
gains derived in ref. 1 for [cj and [g])- The transfer function used in 
equation (19) is expressed in terms of ik in order to maintain the dependence 
of matrix [U^ (eq. (11)) on M and k only. Furthermore, equation (19) 
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represents the overall transfer function and implies that any actuator transfer 
functions should be compensated for to obtain the transfer function matrix. 

On the basis of results obtained in reference 1 , the aerodynamic eigen- 
values are expected to increase if all the gains in matrix [g] are increased 
by a constant factor. Hence, constraints are imposed on the gain values. 

The range chosen for [g] is -100 S G^j i 100. The effects of changing the 
nominal range are investigated in the present report. To reduce to a minimum 
the activation of the control surfaces at zero frequency, the following con- 
straints are imposed on [c] : 


Cii = C 21 = C-|2 - 0 


-2 ^ C22 = 2 

The gain C 22 is varied within a constrained range because of its large 
effects on the aerodynamic eigenvalues, as was shown in reference 1. 

Results for tralllng-edge control .- The trailing-edge control surface was 
first optimized at M = 0.9 on the basis of the curves for ^min against 1/k 
by following the procedure described earlier in this work (with zero weight 
for X niay ) . The optimum gains obtained for [c] and £g] are 


0 0 

0 - 1.76 


[G] = 


The resulting ^min ^max variations with 1/k are shown in figure 3- 

Although the values of ^min improved over those with no control system 

(fig. 2), these improvements are accompanied by a substantial reduction in 
Xmax over most of the 1/k range, with the exception of Improved values at 
the low end of the 1/k range. (Typically, at lower Mach numbers the values 
of ^min less than those shown in fig. 3.) Optimum values of [g] are 

small compared with the wide range of variation permitted the_ Gj^j gains; 
that is, -100 S Gj^j ^ 100. It may be argued that although Xjjj^^ seems to be 
bounded and its values cannot be improved beyond a certain_ maximum level, 

Xmax might be increased without substantially degrading Xjjjj^j^ in accordance 
with_ the relaxed energy requirements. The incorporation of the weighted area 
of" ^max against 1/k into the target function (in addition to the ^min) 
shows that both G 21 and G 22 are positive and increase monotonically with 
the increase in weight associated with ^max* Consequently, an appropriate 
weight could be assigned to ^max which causes the larger term (G 21 ) to reach 


0 0 
4.6 1.53_ 
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the constraint. However, a simpler method is to assign G2^ the value at the 
constraint and optimize for ^rnin only. By following this procedure the opti- 
mized gains for [c] and [g] obtained at M = 0.9 are 

0 0 " 

0 0.46 

’ 0 0 ‘ 

100 79 


The variations of ^max with 1/k are shown in figure 4. By com- 

parison with figure 2 (no controls), the quantity ^min shows improvement over 
most of the 1/k range. In comparison with figure 3, the values of 
figure 4 are mostly negative and degraded. The values of ~^xnax figure 4, 
however, are greatly improved over those in figures 2 and 3 especially^ at low 
values of 1/k. The values of much smaller than those of 

accordance with the objectives set by the relaxed energy approach. 

Effect of variation of optimization range of The results in fig- 

ures 3 and 4 relate to a Gj^j range of 100. Computations made for smaller 
ranges of Gj^j show that the value of C22 becomes increasingly important as 
the Gj^j range is reduced. At the limit, when all the Gj^j terms are set to 
zero, the C22 term provides the only means for improving the eigenvalues. 
(When the Gj^j terms are large, C22 t>e ignored.) Hence, in an attempt 

to reach optimized control values that do not exhibit sensitivities to the 
optimization range of except for a constant factor common to all the 

Gj^j terms, C22 was assigned a value of C22 = -1.86. This value was 
obtained by optimizing the control gains with zero gains. Reoptimizing 

with C22 ° -1.86 and G2^ = 100 resulted in the following transfer function: 




[C] 



0 

- 1.86 


[G] 


0 

100 


0 

80 


Regardless of the value of the G^j constraint, the ratio of G21/G22 remains 
essentially constant. Therefore, the following form of the optimized transfer 
function (eq. (19)) is suggested: 


15 



[t] 


0 



0 


-1.86 J 


+ 


ax(ik) 


0 

100 


0 

80 _ 


( 20 ) 


The free parameter a^ permits the variation of while keeping the ratio 

of G2 i/® 22 ® constant. Figure 5 shows a comparison made at M = 0.9 between 
the curves of X ^-j n against 1/k obtained by optimizing the gains of [T] 
for different maximum values of G21 (fig. 5(a)) and the curves obtained by 
using equation (20) with values of ax chosen to obtain the same maximum 
values for G21 (fig. 5(b)). The values of X ^i n for’ tbe low 1/k range are 
not shown in figure 5 in order to scale the ordinate to yield a good resolution 
among the different curves. (A similar representation, which will be adopted 
in many of the figures, has no reflection on the k range for which the opti- 
mization has been performed.) The differences between the curves in 

figures 5(a) and 5(b) are not significant in view of the large similar values 
of X niav in figures 6(a) and 6(b). These results justify the form of the 
transfer function in equation (20). 


Mach number effects: The effect of changing the flight Mach number on 

Xmin is shown in figure 7 for G21 = 100 and C22 = -1-86. Figure 7(a) is 
obtained by optimizing the gains at the different Mach numbers, whereas fig- 
ure 7(b) shows the values obtained by using equation (20) with a^ = 1 • Once 
again the differences between the X „,^- y, curves in figure 7 are not significant 
in view of the similar and large values of X ^a y shown in figure 8. 

Results for leading -edge — trailing-edge control system.- An identical 
procedure to the one discussed in the trailing-edge system has been adopted for 
the leading-edge — trailing-edge control system in an attempt to optimize X mi n 
and take account of X may . Here again C22 is assigned a value of -1.86 and 
^21, max i® assigned the nominal value of 100. The following optimized gains 
are obtained for the leading-edge — trailing-edge system: 


[C] 



0 

- 1 . 86 _ 


[G] 


-100 

100 


100 ~ 

80 _ 


The variation of X^^jj and X max with 1/k, at M = 0.9, is shown in 

figure 9. The large improvements in ^min> ss compared with the trailing-edge 
system (fig^ 7), can be observed in figure 9(a) for all values of 1/k. The 
values of Xm^x i*^ figure 9(b) are essentially the same as those obtained for 
the trailing-edge system (fig. 8). This situation implies that the optimum 
leading-edge control gains essentially affect with negligible effects 

on Xmay. On the other hand, optimum trailing-edge control gains mainly influ- 
ence X may and have little effect on X mi ^ . 
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Effect of variation of the optimization range of Optimum control 

laws were derived for different values of optimization ranges by using a pro- 
cedure identical to that used for the trailing-edge control system. The 
results obtained for X mi n are shown in figure 10(a) and are compared in fig- 
ure 10(b) with the values of X mi n obtained by using the suggested values 
for [c] and [g] . (The suggested values of [g] are scaled to span the 
appropriate G^j range.) A similar comparison is made in figure 11 regarding 
the variation of Xmav The agreement obtained in figures 10 and 11_ together 
with the localization of the leading-edge effects to the values of X^^n lead 
to the following more general formulation of the control gains; 



0 0 




0 

0 

0 

l_L 

= 


+ ik 





0 -1.86 


a-p 


100 80 


where sl and a^ are positive free parameters. The value of Sl essen- 
tially determines the levels of X^jl^ whereas the value of a'p essentially 

determines the levels of X ^a v • If the value of aL = 0 is substituted in 

equation (21), the resulting matrix [t] is identical to the one pertaining 
to the trailing-edge control, as presented in equation (20). Hence, equa- 
tion (21), together with the parameters Sl and aip, is applicable to both 
trailing-edge and leading-edge — trailing-edge control systems. 

Mach number effects: The effects of Mach number on the variation of X ^j p 

with 1/k for G 21 = 100 and C 22 = -1.86 are shown in figure 12. Fig- 
ure 12(a) is obtained by optimizing the gains at the different Mach numbers, 

whereas figure 12(b) shows the values of Xjj,^^ as a function of 1/k as 
obtained by using equation (21) with aj^ = a-p = 1 • As can be seen, the fig- 
ures are almost identical. Similar results are obtained for the accompanying 
variation of with 1/k shown in figures 13(a) and 13(b). Figures 14 

and 15 are identical to figures 12 and 13 except that they show the variation 
of the different X values at the low end of the 1/k range. A deterioration 
in the values of Xjjjj^jj and X ,^ay occurs at low 1/k values as the Mach num- 
bers are reduced and especially when M = 0. However, since the results per- 
taining to M = 0 are meaningless for practical applications and since the 
reduced values of at the high 1/k range are shown in the following 

section to be relatively less important, no attempt is made to improve the 
Xmin behavior at low 1/k values at the expense of some degradation at the 
higher range of 1/k. 

Effect of structural damping .- To understand the meaning of a specific 
value of an aerodynamic eigenvalue, it is convenient to express the structural 
damping in terms of an equivalent aerodynamic eigenvalue X^j (the full devel- 
opment of the expression for X^j is given in appendix B) . The structural 
damping force acting on the ith mode is given by equation (B1) as 

2 

Fd,i = iSi‘^n,ibRmiqi (22) 
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where is the structural damping coefficient, Wn,i is the natural fre- 

quency, m^ is the generalized mass, Qi is the nondimensional generalized 
modal response, and bp is a reference wing semi chord. The subscript i 
relates the parameters to the ith generalized mode. It can be shown (see 
appendix B for details) that the work W done per cycle by the force Fd,i 
is given by equation (B7) as 


j 2 2 2 

W = — Tf^pbg^V^SQXjj^ j^qQ^i = , i^R^™i^o , i 


where bp is the semichord at the control-surface midspan region and Sq is 
the control-surface span. This equation yields (eq. (BIO)) 


X 


d,i 



(23) 


where is defined by equation (B11) as 


*^n, i^c 


^n,c,i 


(24) 


For a two-dimensional strip, s n Sg, bp = bg, and equation (23) reduces to 


^d,i = , c , i^R , i 


(25) 


Equation (25) shows that a fixed value of structural damping converts to an 
increasingly larger equivalent aerodynamic eigenvalue as the reduced fre- 
quency k„ c i is increased. Therefore, a value of X^iji = -1, for example, 
at kjj Q j^’ = ’o.1 is equivalent to a structural damping coefficient 1225 times 
larger than the one obtained if the same value of X^j^n = -1 is obtained at 




3.5 


^1225 



This example illustrates the relative importance 


of the aerodynamic eigenvalues in different regions of the k range and 
stresses the importance of the low k range over the higher k range. 


Summarizing remarks regarding damping type transfer function.- The 
relaxed energy approach, when used in conjunction with the damping type trans- 
fer function, yields fixed control gain ratios for both the trailing-edge and 
the leading-edge — trailing-edge control systems, which are effective over a 
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wide range of reduced frequency k and of control gain factors a^ and a-p. 

It is shown that although the trailing-edge control system shows improvements 
in the values, these improvements are bounded and cannot be improved 

beyond a certain maximum value. The values for Xju^x? however, are shown to 
be unbounded and increase continuously with a-p. The leading-edge — trailing- 
edge system shows unbounded improvements over a wide range of k values for 
both Xjjjj^j^ and X^ax* gain ap, essentially affects the Xju^^ values 

with only little effect on X^ax? whereas the gain a«p essentially affects the 
values of with little effect on X^pj^. This unbounded behavior is typi- 

cal of a dynamic system to which damping is added. For this reason the [g] 
matrix is regarded essentially as a damping matrix introduced through the acti- 
vation of the control surfaces. 

Two major criticisms can be made regarding the damping type transfer 
function: 

(1) The damping type transfer function does not permit the introduction 
of damping into the system over a limited range of k values. Instead, damp- 
ing is introduced over an uncontrolled range of k values. Hence, a penalty 
may be paid in terms of unnecessary control-surface activity. 

(2) The term ik[o] that appears in the damping type transfer function 
(eq. (19)) increases indefinitely with k. As a result, stability problems 
may develop at very high k values beyond the range investigated in this work. 

In addition to these criticisms, the relaxed energy approach does not rule 
out the use of stiffness terms in order^to change the response of the system to 
increase the response associated with X^ax- damping type transfer func- 

tion permits the introduction of aerodynamic stiffness terms through the matrix 
[c] . This situation is inconvenient and very limiting since the matrix 
has a direct effect on the static behavior of the system. 

These points lead to the formulation of the localized damping type trans- 
fer function which is presented later. 


Localized Damping Type Transfer Function 

The following transfer function is a direct outgrowth of the criticism of 
the damping type transfer function: 



(ik)2 

^ 

(ik)2 + 2^kn(ik) + k„^ 



( 26 ) 


where both 5 and k^^ are constants. Direct comparison between equa- 
tions (26) and (18) shows that the term ik in the damping type transfer func- 
tion has been replaced by the term (ik)^/[[(ik)^ + 2?kj^(ik) + k^^^J . It is 
appropriate, therefore, to investigate the properties of this term and deter- 
mine how it compares with the ik term. 
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Variation with k of localized damping type transfer function.- The vari- 
ation with k of the localized damping type transfer function (L.D.T.T.F.) is 
essentially embedded in the term (see eq. (26)) 


(ik)2 

R = 

(ik)2 + 2?kjj(ik) + k„2 


(27) 


Defining X as 


k 

X = — 


and substituting in equation (27) yields 


( 28 ) 


-X2 

R = — 

-X2 + 2iCX + 1 


(29) 


The denominator in equation (27) is identical in form to a transfer function 
representing a second order system with damping C and natural frequency kj^. 
Figure 16 shows the variation of the real and imaginary parts of R as a 
function of X and C. It can be seen that R = 0 when k = 0 and R h 1 
when k ^ The imaginary part of R is always positive. For small values 
of very large values of Rj (that is, the imaginary part of R) are 
obtained, with a maximum value obtained around k = k^^. In addition, the vari- 
ation of Rj with k is very rapid for small values of As the value 

of C is increased, Rj becomes smaller and its maximum values occur for 
k > kjQ. In addition, the increased value of C reduces the variation of Rj 
with k. If Rj is associated with positive damping, then positive stiff- 
ness values are associated with positive values of Rp (that is, the real 
part of R) and negative stiffness values are associated with negative values 
of Rp. (Strictly speaking, some symmetrical aerodynamic stiffness terms may 
be introduced through formal multiplication of the aerodynamic forces and [t] . 
However, these terms appear to be small.) Therefore, the quadrature term ik 
can be locally simulated by R over a range of frequencies determined by 
both C and kj^. In addition, stiffness terms are introduced as R varies 
with k. These stiffness terms clearly vanish when k = 0 and therefore do 
not affect the static behavior of the system. A transfer function of the 
type ik/(ik + a), where a is a real number, could also simulate the quadra- 
ture term ik, however, with reduced effectiveness and reduced control regard- 
ing "localization" of damping. 

Remarks on localized damping type transfer function .- Because of the 
results obtained for the damping type transfer function and the similarities 
between the two transfer functions, the following values have been adopted for 
the localized damping type transfer function: 
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(1) C22 is assigned the value of -1.86 (derived for zero gains and 

therefore also applicable to the localized damping type transfer function) to 
maintain the effectiveness of the system for small, values of matrix [g] . (All 
other Cij values are maintained at zero.) 

(2) A nominal optimization range for the Gj^j elements is chosen as 

-4 ^ Gij S 4. The smaller_ range for G^j was chosen to yield approximately 
the same peak values for X as obtained for the damping type transfer function 
at the low and intermediate ranges of k values. 

Results for trailing -edge control system .- The trailing-edge control sys- 
tem is first optimized at M = 0.9, Gj^j being constrained within the nominal 
range, and with kj^ = 0.2 and 5 = 0.5. Variation of these values is subse- 
quently investigated. 

Initial optimum values obtained for the matrix [g] showed behavior 
identical_ to that of the damping type transfer function; that is, a deteriora- 
tion in X mav was obtained, and the resulting Gj^j gains were very small in 
comparison with the span of the nominal Gj_j range. By following an identical 
type of investigation as in the damping type transfer function case, an iden- 
tical conclusion was drawn, that is, that optimization must be performed after 
assigning G2i the maximum positive value of its range. The results for the 
nominal range yield the following gains: 


0 

- 1.86 


0 

2.8 


Figure_17 shows the resulting variation of Xj^j^n and Xj^^x with 
Both Xjnin ^max positive throughout the 1/k range and 

shows "peaking" behavior_at kjj (that is, 1/k = 1/0.2 = 5). Figure 18 
shows the variation of Xm^^ and X ^av with 1/k at different Mach num- 
bers_ obtained by using these gains. Figure l8(a) shows a slight deterioration 
in Xmj^n as the Mach number is reduced and shows no signs of peaking at kjj. 
Figure l8(b) shows that the peaking characteristic of X ma x _exists at all 
values of M. Figure 19 shows the variation of and X ^a v with 1/k at 

M = 0.9 obtained by using these gains and scaling the [g] matrix by various 
constant factors. The results show increased X ^i^y values throughout the 
range for most of the scaled gains of the [g] matrix. A small variation in 
Xmin seen in figure 19(a), especially for the case where these values 

of [g] are all multiplied by 2. Figure 20 shows the variation of and 

Xjjiax with C . The effect of varying ? is seen to be similar to the effect 
of scaling [g] (fig. 19), with the exception that the X ^ gy peaks become 
narrower as ^ is reduced and the dips which follow the peaks become more 
accentuated. Figure 21 shows the variation of and X ^ay with 1/k 


1/k. 

^max 
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are 


at M = 0.9 for different values of kjj. The large improvements in X ^av 
immediately apparent over a wide range of k values. The values of ^niin> 
however , show a negative region for the value of k„ = 0 . 1 . Comparison with 
the curve for the unactivated system (fig. 2) shows that around k = 0.3 

= 3.5 where X^jj^jj is smaller, the degradation in Xj^^^ is very small. 

The following generalized control law is therefore suggested for the 
trailing-edge system: 



0 1 ay(ik)2 

+ ■ 

•1.86J (ik)2 + 2i?knk + kn^ 



(30) 


where a^, C, and k„ are left as free parameters which determine the amount 
of damping introduced, the distribution as a function of k, and the location 
of the peak along the k-axis. 

Results for le ading-ed ge — trailing-edge control system.- The following 
optimum gains are obtained for the leading-edge — trailing-edge control system 
at M = 0.9 by using the nominal range for the gains (with k„ = 0.2 

and C = 0.5): 


[c] 



0 

- 1.86 


[G] 



4 

2.8 


The resulting variation of Xj^^jj and X ^ay with 1/k is shown in figure 22. 
The variation of X ^^y with 1/k is identical to that obtained for the 
trailing-edge control system (fig. 17) and implies that the leading-edge con- 
trol gains have little or no effect on X^j^x- It has been shown that the 
trailing-edge control gains associated with [g] have only a small effect on 
Xmin- Therefore, the large values of X^in shown in figure 22 (a2 are essen- 
tially due to the leading-edge control gains. The variation of and X ma v 

with Mach number, using these control gains, is shown in figure 23. Once 
again, X ^g y remains essentially the same as that for the trailing-edge sys- 
tem whereas ^min is greatly improved and shows large peaks at k = k^ = 0.2 
(or 1 /k = 5 ) . 

The variation of the aerodynamic eigenvalues with the_ scaling factor of 
the matrix [g] is shown in figure 24. Both Xjjj^n and X^j^x increase with 
the scaling factor over most of the range; this result implies that both X^j^n 
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and Xmax unbounded in the leading-edge — trailing-edge system. The 

unbounded characteristic of the eigenvalues is observed once again in figure 25 
where the value of C is varied. The peaks of X^^x figure 25 are essen- 
tially located at k = kjj whereas Xj^j^^j shows a peak at this location only 
for small values of C . As the value of ? increased, the peaks move toward 
lower values of 1/k. A study_of the peak values in conjunction with R 
(fig. 16) indicates that the X ^av peaks occur at Rr = 0, or k = k^. On the 
other hand, the Xm^^ peaks occur at the value of k which yields the largest 
response in quadrature Rj. Figure 16 shows that as C increases, Rj ma y is 
located along curves for which k > kjj. Therefore, if it is desired to maxi- 
mize simultaneously X^^jj and X^j^x a given value of k, the leading-edge 

control should be given a similar transfer function but with a smaller value 
of kfj. The effect of kjj on the aerodynamic eigenvalues is shown in fig- 
ure 26. The results conform with the foregoing discussion. On the basis of 
these results, the following generalized control law is suggested; 


{ 





aL(ik)2 

0 

\ 

To 0 

0 -1.86 

+ 

(ik)2 + 2CLkn,L(il«) + ^n,L 
0 

ai-(ik)2 

-4 4l 

4 2.8 1 


(ik)2 + 2^Tkn,T(iX) + k^^T 



(31) 

where the subscript L refers to the leading-edge control and the subscript T 
refers to the trailing-edge controls. The six free parameters allow the ampli- 
tude of the different X distributions, the width of the distributions, and 
their peak locations to be selected. If a^ = 0, equation (31) reduces to the 
trailing-edge control law given by equation (30). 

Summarizing remarks regarding localized damping type transfer f u nction . - 
The localized damping type transfer function is shown to have basically the 
same characteristics as the damping type transfer function regarding its 
effects on the values of ^max».^^^ within a localized range of 

k values. The deteriorating effects on shown by the damping type 

transfer function at high k values have either been eliminated or greatly 
reduced. Further improvement in ^min high k values can be obtained by 

means of the free parameters left in the generalized control law (eq. (31)). 

In addition, it can be shown that at sufficiently high_ values of k, the aero- 
dynamic eigenvalues are proportional to k, that is, \ ^ ky whereas the 
equivalent aerodynamic eigenvalue for the structural damping behaves 

as X^^ oc k^. Therefore, at sufficiently high frequencies, X^ is dominant 
and the system is always stable (not the case for the D.T.T.F. where X ^ k^) . 
Since the values of highest value of k considered ( = 3.5) were 

seldom negative and assumed relatively small values, no special reference to 
this region was made during the discussion of the results. Finally, the free 
parameters left in the generalized control law, together with the stiffness 
properties exhibited by the localized damping type transfer function, make it 
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possible to tailor fit the results obtained herein to any specific system with- 
out exposing it to sensitivity hazards. 


Sensitivity Test for the Gains 

Sensitivity tests were conducted on the and Gj^j gains for the two 

transfer functions and the two control systems investigated in this work. All 
the Gij terms were found to be important (either for Vin ^max)» 

and no sudden deteriorations were observed with small to moderate parametric 
changes. Figure 27 is an illustrative example in which the value of G22> 
using a trailing-edge system with localized damping type transfer function, 
is varied at M = 0.9. Both X values are relatively insensitive to small 
changes in G22- Finally, the gains obtained by the relaxed energy approach 
are different only in magnitude but not in sign as compared with those obtained 
in reference 1 . 


SIMPLIFIED CONTROL LAWS 

The suggested generalized form of the damping type transfer function is 
given by equation (21) as 


~0 0 " 


^L 


-100 100' 


+ ik 




0 -1.86 


a.*p 


100 80 


Equation (21) can be simplified if a constant reference value of velocity Vp 
is assumed in the expression for k. In this latter case, the only variable 
is w and therefore [t] can be written in the form 


[T] 



"o 0 

(A) 

aL 


“-100 100' 

= 


+ i — 





0 -1.86 

0)R 

a'p 


100 80 


(32) 


Vr 

where wr = — is a reference frequency. Since aL and a^ are free param- 
b 

eters, it is convenient to reduce the scaling of the element by a factor 

of 25 to make the maximum Gj^j values the same as for the case of the local- 
ized damping type transfer function. Hence, equation (32) (scaled) yields the 
following standardized form of control law for the damping type transfer 
function; 



(33) 
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(Note that for a given case the values of sl and a^ in equation (33) differ 
from those in eq. (32).) Similarly, if 


0)nb 



the localized damping type transfer function can be simplified to yield the 
following form of standardized control law (see eq. (3D): 


{ 


e 

6 


} 


/ 

aL(iw)^ 

0 

\ 

[0 0 


(ito)2 + 2^L“n,LDw) + “n,L „ 

-4 4I 


+ 

aT(ito)2 


CO 

I 

0 


0 ; ^ 2~ 

4 s.sji 


(ico)^ + ^ «p( i(o) + 

/ 



(34) 


The similarity between the gains appearing in equation (33) and those appearing 
in equation (34) are now immediately apparent. 


If aj__ and a^ are sufficiently large to make the effect of C 22 negli- 
gible, then equation (34) shows that the ratio between h/b gains and the 


a gains is 1:0.7 



Since these gains are associated with a point at 0.6b 


from the leading edge (30-percent chord), identical activation will be obtained 
by a linear sensor at 1.3b (=0.6 + 0.7) driving the trailing-edge control sys- 
tem. This point is at the 65-percent chord measured from the leading edge (the 
damping type transfer function yields a point located at 70-percent chord). 

For either transfer function the leading-edge control can be driven by a single 
sensor located at -0.4b (=0.6 - 1), or 20-percent chord ahead of the leading 
edge. This location is clearly not practical and, therefore, two accelerom- 
eters should be employed to drive the leading-edge control . 


These generalized control laws essentially determine the relation between 
the gains associated with the h/b and a for each control surface irrespec- 
tive of whether one or more control surfaces are actuated on the same strip. 
This statement implies that the relationship between the gains will be appli- 
cable to three-dimensional flows, if the ratios between the lift and the moment 
produced by the rotation of the respective control surface are the same as in 
the two-dimensional flow case. 


CONCLUDING REMARKS 

Control laws using realizable transfer functions are derived which per- 
mit the relaxation of the stability requirements of the aerodynamic energy 
approach. The resulting aerodynamic eigenvalues indicate that both the 
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trailing-edge and the leading-edge — trailing-edge control systems can be made 
more effective. As a result, flutter suppression and gust-alleviation problems 
can successfully be tackled by the aerodynamic energy approach with a trailing- 
edge control system only and without encountering sensitivity to changing 
flight conditions. In the original aerodynamic energy approach, such an 
ability could be handled by a leading-edge — trailing-edge system only. 

Free parameters have been identified which control the amount of the aero- 
dynamic damping introduced by the activated system. Furthermore, one of the 
control laws also permits the selection of a reduced frequency range within 
which damping is to be increased and provides stiffness terms for controlling 
the response of the system. The values of these free parameters can be 
tailored to specific applications. The results have been reduced to a simple 
form and are applicable to a very wide class of aircraft operating within the 
subsonic Mach number range. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 7, 1977 
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APPENDIX A 


THE ENERGY ANALYSIS 


Let the n equations 

{f} = + upb^sCAp + iAj)J {q} + [e] {q} (A1) 

represent the equations of motion of n structural modes with r activated 
controls, where at flutter 

{f} = 0 


and 

0) 

[B] 

M 

[E] 

P 

s 

b 



frequency of oscillation 
mass matrix 

real and imaginary parts of aerodynamic matrix, respectively 
stiffness matrix 
fluid density 
reference length 
reference semichord length 


All matrices in equation (A1) are of the order n x (n + r) (n structural modes 
and r active controls). The response vector {q} in terms of n structural 
responses and r control deflections is 


{q} 


qc 


(A2) 


Equation (A1) can therefore be written as 



APPENDIX A 


Assume a control law of the form 


{qo! = M to) 


(A4) 


where [t] is a r x n matrix representing the transfer functions of the con- 
trol law. By substituting equation (A4) into equation (A3), the following 
equation is obtained: 


{F} = (-«i2[[B3] . [Bj[T] . Fpb‘«s([A„,3] . [A„,J[I] 

[‘i,s] * M - [Ec]M){p) 


+ 1 


(A5) 


where the products [BqJ [t] , [Ar,c3 1X1 > [*I,c] [X > E^c] [X 

n X n matrices. The vector {q} can be written as 


{q} = = {qR + iqilei^t 


(A6) 


The real part of equation (A5) is essentially the same as the imaginary 
part, except for the initial conditions at t = 0. In harmonic motion, where 
no transients are treated, only the real part of equation (A1) need be con- 
sidered. Hence, the forces which the system exerts on its surroundings are 
given by 


{f} = 


M + [Bc][X ^ ^Pb"*s([AR^3] ^ [AR,e][X i[Al.s] + i[Al,c][X) 


^ fc] ^ XcjrX {qj®^"^^ [Bs] ^ [Bc][T»] 


^PB'*s([Ar,s] + [Ar,c][t»] - i[Ai,s] - i[Ai,c][T*]) 


+ [Eg] + [Ec][T*]j{qo*}e-i“t (A7) 

where the asterisk denotes the conjugate vector. The vector {f} is clearly 
a real vector. 
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APPENDIX A 


The velocity vector {q} can be obtained by differentiating equation 
to obtain 


{q} = iO){qo}ei‘^t 

The real part of the velocity vector R(q) is given by 

R{q} = -1 iwiqoleif^t _ 1 io){qo»}e-i“t 


Hence, the rate at which the system does work on its surroundings is given 
by K[qJ{p}. Equations (A?) and (A9) yield 


R[_qJ{F} = [Bs] [Bc][t] + TTpb'^s ([Ar , 3 ] + [Ar,c] M 


- i[Ai,s] - i[Al,c][T])] - [Es] - [Ec][T]]{qo}e2i^t 

- " ^Pt^^3([AR,3] . [Ar,J[t] 


+ 1 


‘[*I,s] ■> i[''l,o][T])] ^ Ps] - M[T])(qol 


* * MM * "Pb‘'s([AB,s] * [Ab.JM 

- - i[Al,o][T«])] - [Es] * [Eo] [T*]) {q/} 

- rMj(-“^[M * MM » l-Pb"q([A„,J . [Ar.JM 

- l[Ai.J - l[Al,o][T*])] . [E 3 ] - [Ec][T*]){q<,»le-21“b 


(A 6 ) 

(A 8 ) 

(A9) 


(A10) 
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Hence, the work P done per cycle by the system on its surroundings can be 
found by integrating equation (A10) between t = 0 and t = 2ir/o). Thus, 

P = * [Bo]M * >tPb‘'s([AR,J + [Ar,o][T*] 

- ^[* 1 , 3 ] - i[Al.o]&»])] * [E 3 ] * [Eo][T»]jko*l 

- ^LPo*J(-“2[[Bs] * [Bc][T] * 1-Pb“3([A„,3] * [Ar,J[t] 

+ i[*I,s] + i[Ai,o]W) * [Es] t [Eo][T]j{qo) 

This equation is a scalar equation; therefore, the first expression on the 
right-hand side of this equation can be transposed ^ and [Eg] being 

assumed to be symmetric^ 

0 = "-^^'KJ[-([Ai„] [Al.a]" * [Al,o][T] * M " [Ai , o] " ) 

* l([»».s] - [Ah,s]" * [Ar,o][t] - [T*]"[Al,o]’')]{qol 

- |K*J [i([Bc] ra - M ^[Ec] ^)] {Po} (A1 1 ) 

The matrix [Eq] is assumed to be zero since, in general, no elastic couplings 
exist between the structural modes and the control deflections. Hence, equa- 
tion (All) reduces to 
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TT^pb^w^s 

p = ^ — L^o*J 


■([«I,3] + [Ol.s]’' + [Al.o][T] + [T*]^[»I,o]^) 


,/[Bo][T] - [T-]"[Bo]" ^ 

^ pR,s] - [Ar,s] 


Tfpb^S 


[Ar,o][t] - [T*]’'[Ar,oJ"' 


{Qo} 


(A12) 


Note that the matrices within the square brackets form a Hermitian matrix. 

Now determine the eigenvalues and eigenvectors of the following Hermitian 
matrix [u] extracted from equation (A12): 


x{n} = [uj{n} 


(A13) 


where 


[u] = 


-([Al,s] - [A1.3]" - [Al,o][T] * pAj^LAl.e]^) 


+ + [Ar,s] - [Ar.s] 

TTpb^S 


- [AR,o][T] - [T*]^[Ar,o]^' 


(A14) 


and represent the vector {qol terms of the eigenvectors of equation (A13), 
that is, 


{Qo} = [Qr + iQ±]{CR + i?l} 


(A15) 


The columns of the square matrices [Qp] and [Qj] are the real and imaginary 
parts of the eigenvectors of equation (A13). The vector + i^j} contains 

the generalized modal coordinates of the transformation defined by equa- 
tion (AI 5 ). 


All the eigenvalue and eigenvector solutions of equation (A13) can be 
expressed in the form [Qr + iQi] |^Xj = [u] [Qr + iQj] • Postmultiplying this 
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equation by {Cr + i^j} and premultiplying it by [Qo*J 
[Cr - iCiJ [Qr"^ - iQi'^] yields 

L^R - [Qr'^ - [Qr + iQl] K{Cr + i?l} 


= L^R - [Qr'^ - 


-([Ai,s] * [*I,s]^ + [*I.o][T] 


T /[®o 1 w - [T*]^r®oi'' 

- * [Ar.J - [ar.=]" 

\ irpb^s 


- [AR,o][T] - [T*]"[A„,,]' 


[Qr + iQi]{CR + i?i} 


(A16) 


where [xj is the eigenvalue matrix, that is, 


Xl 

^2 



(A17) 


The right-hand side of equation (A16) is identical to the right-hand side of 
equation (A12) except for the iT2pb^w^s/2 factor. Hence, the left-hand sides 
of these two equations can be equated to obtain 


TT^pb^W^s P 

P = [Cr - iCiJ [Qr^ - iQi'^J [Qr + iQl]{CR + i?il 


The matrix [Qr"^ - iQj'^J represents the modal matrix of the complex left-hand 
normalized eigenvectors and therefore 


[QrT - iQjT] [Qr + iQj] = [i] 


(A18) 
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where [f] 


is the unit matrix. Hence, 


P reduces to the form 


P = - i5iJPl]{Cb * i5i} 

P = * LfiJNfe}) 


The energy input per cycle into the surroundings has thus been reduced to 
quadratic form in terms of the generalized energy modal coordinates. 

Note that if all control surfaces are mass balanced, the matrix [BqJ 
vanishes and [jjJ reduces to 


[D] = [-(Pl,s] - Pl,s]^ * [Al,o][T] - 

* i([Ar,s] - [Ar,s]" - [Ar.o]W - M"[Ar.c]")] 

when the control system is not activated [t] = 0 and 

M = [-([Al,s] * [Ai,s]") - i([AR,s] - [Ar.s]")] 


(A19) 


(A20) 


(A21) 


(A22) 
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CONVERSION OF STRUCTURAL DAMPING INTO EQUIVALENT 
AERODYNAMIC EIGENVALUES 

Let the structural damping force F(j ^ acting on the ith structural 
mode be expressed by 

Fd,i = iSi“n,it5Rmiqi (B1) 


where 


Si 

“n,i 

“i 

t>R 

Qi 

The work 


structural damping coefficient of ith mode 
natural frequency of ith mode 
generalized mass of ith mode 
reference wing semichord 
nondimensional generalized modal response 

W done by the damping force per cycle of oscillation is given by 


W = 


J -» 2 ir/o) 
0 


P’s 


dqi 

dt 

dt 


(B2) 


Substitution of equation (B1) into equation (B2) yields 


2 o dq 

W = SiWn,ibR^“i ^'5 TT 

oq dt 

Let 


(B3) 




(BH) 


<3i 

iqj^ = — = qoji cos cot 


(B5) 
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! 

■i 


Therefore, 

2 2 

W = TigiWnjibR^miqo^i 


(B6) 


Let the work done per cycle be defined in a form equivalent to the newly 
defined aerodynamic eigenvalues, that is, 



^ 2 2 
Tt2pb^2v2s^Xd,iqo,i = iTgiiOn , ibR^miqo , i 


(B7) 


where is the semichord at the control-surface midspan region, Sq is the 

control-surface span, and ^ is the equivalent value of the structural 
damping in terms of the aerodynamic type eigenvalues. Equation (B7) yields 


2giWn , i*^R^®i 

Trpbc2v2sQ 


(B8) 


Define the mass ratio Mr^i as 


m^ 


Mr,i = 


upbR^s 


(B9) 


Substitution of equation (B9) into equation (B8) yields 


X 


d,i 




(BIO) 


where is defined by 


*^n , c , i 


^n , i^c 


V 


(BID 


35 





APPENDIX B 


For a two-dimensional strip s = Sq and bp = b^. Equation (BIO) reduces in 
this case to 


^d,i = ^ Q ^ ^ i 


(B12) 
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Figure 2.- Variation of X with 1/k at various Maeh numbers for a wing strip 

with no control surfaces. 
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Figure 3.- Optimum variation of X with 
(with gains optimized on the basis of 
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(b) ^max* 

1/k for a trailing-edge control system at M = 0.9 
Xmin only) using damping type transfer function. 
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(a) Optimized gains. (b) Fixed gains. 

Figure 5.- Variation of with 1/k using optimized gains for different [Cijl ranges 

compared with a similar variation using fixed gains scaled by aip. Trailing-edge control 
system at M = 0.9 and using damping type transfer function. 




(a) Optimized gains. (b) Scaled gains. 

Figure 6.- Variation of X mav with 1/k using optimized gains for different ICijl ranges 
compared with a similar variation using scaled gains, trailing-edge control system at 
M = 0.9, and damping type transfer function. 
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(a) Optimized gains. 


(b) Fixed gains. 


Figure 7.- Variation of X ^-j n with 1/k using optimized gains at different Mach numbers 
compared with a similar variation using fixed gains, trailing-edge control system, and 
damping type transfer function. 
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(a) Optimized gains. (b) Fixed gains. 

Figure 8.- Variation of with 1/k using optimized gains at different Mach numbers 

compared with a similar variation using fixed gains, trailing-edge control system, and 
damping type transfer function. 
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Figure 9.- Optimum variation of X with 1/k for a leading-edge — trailing-edge control system 
at M = 0.9 (with G 2-1 = 100) using damping type transfer function. 
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(a) Optimized gains. (b) Scaled gains. 

Figure 10.- Variation of Xmin with 1/k using optimized gains for different |Gijl ranges 
compared with a similar variation using scaled l^ijl gains, leading-edge — trailing-edge 
control system at M = 0.9, and damping type transfer function. 
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(a) Optimized gains. (b) Scaled gains. 

Figure 11.- Variation of with 1/k using optimized gains for different |Giil ranges 

compared with a similar variation using scaled l^ijl gains, leading-edge — trailing-edge 
control system at M = 0.9, and damping type transfer function. 
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(a) Optimized gains. (b) Fixed gains. 

Figure 12.- Variation of ^min with 1/k using optimized gains at different Mach numbers 
compared with a similar variation using fixed gains, leading-edge — trailing-edge control 
system, and damping type transfer function. 
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(a) Optimized gains. (b) Fixed gains. 

Figure 13.- Variation of ^max with 1/k using optimized gains at different Mach numbers compared 
with a similar variation using fixed gains, leading-edge — trailing-edge control system, and 
damping type transfer function. 
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(a) Optimized gains. (b) Fixed gains. 

Figure 14.- Variation of with 1/k using optimized gains at different Mach numbers compared 

with a similar variation using fixed gains, leading-edge — trailing-edge control system, and 
damping type transfer function. 
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(a) Optimized gains. (b) Fixed gains. 

Figure 15.- Variation of with 1/k using optimized gains at different Mach numbers compared 

with a similar variation using fixed gains, leading-edge — trailing-edge control system, and 
damping type transfer function. 
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(a) ^max* 

Figure 17 •- Optimum variation of X with 1/k for a trailing-edge control system 
at M = 0.9 (with G 21 = 4) and using localized damping type transfer function. 
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Figure 18.- Variation of X with 1/k at different Mach numbers. Trailing-edge control system 

using localized damping type transfer function. 
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Figure 19.- Variation of X with 1/k at M = 0.9 for different scaled ranges of jG^jl 
Trailing-edge control system using localized damping type transfer function. 



(a) ^min* ^max* 

Figure 20.- Variation of X with 1/k at M = 0.9 for different values of 
Trailing-edge control system using localized damping type transfer function. 
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Figure 21.- Variation of X with 1/k at M = 0.9 for different values of 
Trailing-edge control system using localized damping type transfer function. 
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Figure 23.- Variation of X with 1/k at different Mach numbers. Leading-edge — trailing-edge 
control system using localized damping type transfer function. 
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(a) ^min* ^max* 

Figure 24.- Variation of X with 1/k at M = 0.9 for different scaled ranges of |Gij|- 
Leading-edge — trailing-edge control system using localized damping type transfer 
function. 
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Figure 25.- Variation of X with 1/k at M = 0.9 for different values of 
Leading-edge — trailing-edge control system using localized damping type 
transfer function. 
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(a) ^min* ^max* 

Figure 26.- Variation of X with 1/k at M = 0.9 for different values of kjj. 
Leading-edge — trailing-edge control system using localized damping type 
transfer function. 
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Figure 27.- Effects on X of variation of G 22 around its optimum value. Trailing-edge 
control system at M = 0.9 using localized damping type transfer function. 
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